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ABSTRACT 

We analyzed XMM-Newton data of the Perseus cluster out to ~1 Mpc, or approximately half 
the virial radius. Using the flux ratios of Lya lines of H-like Si and S to Ka line of He-like Fe, 
the abundance ratios of Si/Fe and S/Fe of the intracluster medium (ICM) were derived using the 
APEC plasma code v2.0.1. The temperature dependence of the line ratio limits the systematic 
uncertainty in the derived abundance ratio. The Si/Fe and S/Fe in the ICM of the Perseus 
cluster show no radial gradient. The emission- weighted averages of the Si/Fe and S/Fe ratios 
outside the cool core are 0.91 ± 0.08 and 0.93 ± 0.10, respectively, in solar units according to the 
solar abundance table of lLoddersI (l2003h . These ratios indicate that most Fe was synthesized by 
supernovae la. We collected -ftT-band luminosities of galaxies and calculated the ratio of Fe and Si 
mass in the ICM to -ftf-band luminosity, iron-mass-to-light ratio (IMLR) and silicon-mass-to-light 
ratio (SMLR). Within ~1 Mpc, the cumulative IMLR and SMLR increase with radius. Using 
Suzaku data for the northwest and east directions, we also calculated the IMLR out to ~ 1.8 
Mpc, or about the virial radius. We constrained the SMLR out to this radius and discussed the 
slope of the initial mass function of stars in the cluster. Using the cumulative IMLR profile, we 
discuss the past supernova la rate. 



Subject headings: galaxies:clusters:individual(the Perseus cluster)- X-rays:intracluster medium 



1. Introduction 

Clusters of galaxies are the largest gravitation- 
ally bound objects in the Universe. The intra- 
cluster medium (ICM) contains a large amount of 
metals, synthesized by supernovae (SNe) in clus- 
ter galaxies. Thus, the distribution of metals in 
the ICM provides important information on the 
chemical history and evolution of clusters. 

Because metals were synthesized in galaxies, 
the ratios of metal mass in the ICM to the 
total light from galaxies in clusters or groups, 
(i.e., metal-mass-to-light ratios) are key param- 
eters in investigating the chemical evolution 
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of the ICM. Using Ein stein and Gin ga data, 



Arnaud et al.l (|1992D and iTsuru I (|1992f ) derived 
ratios of Fe mass in the ICM to the total light 
from galaxies, which is the iron-mass-to-light ra- 
tio (IMLR). To account for the observed IMLR, 
either the past average rate of SN la was at 
least higher by a factor of 10 than the present 
rate or massive stars in clusters formed with 



a ver y flat initial mass function (Rc nzini et al 
1993). With ASCA observations, the derived 



IMLR within a radius where the ICM density 
falls below 3 xlO -4 cm -3 is nearly constant in 
rich clusters and decreas es toward poorer systems 
(Ma kishima et al.l [200lh ■ With Suzaku, Chan- 



dra, and XMM satellites, the IMLR of several 



measured out to 0.2ri8n-0.5risn ( 


Matsushita et al. 


2007a; Komivama et al.l 2009; 


Sato et al.l 


2007. 


2008L |2009a|b 


. 2010; Rasmussen & Ponman 


2009; 


Sakuma et al. 


2011) and that of the Coma clus- 



ter was derived out to 0.5ri 80 ([Matsushita et al 



1 



2012). Suzaku first measured the Fe abundance 



of the ICM beyond .5riso dFuiita et all I2 OO8: 



Simi onescu et al.ll2011l ). With Suzaku. ISato et al.. 
( 20121) derived the IMLR profile of the Hydra A 
cluster out to rigQ. Here, riso is the radius in 
which matter at 180 times the critical density of 
the universe is contained. In individual clusters, 
the IMLR is lower around the center, indicating 
that Fe in the ICM extends farther than stars. 
Therefore, to derive the total Fe mass in the ICM, 
we need observations out to the virial radius. 

Since Fe is synthesized by both SNe la and 
by core-collapse SNe (hereafter SNecc), we need 
measurements of abundances of various elements 
to constrain contributions from the two types of 
SNe. The ASCA sate llite first studied the S i abun- 
dance in the ICM (iFukazawa et ail Il998l l2000t 



Fukazawa et al 



Finoguenov et al.l 120001 120011) 
il99J) reported that the Si/Fe rati o in the ICM in- 
crease s with ICM temperature, and lFinoguenov et al 



(2000) reported that the Si/Fe ratio increases 
with radius in several clusters. Using Chandra 



data o f groups out to rsoo , Rasmussen fe Ponmanl 



(|2007l ) reported that the SNcc contribution in- 
creases with radius and completely dominates at 
rsoo- XMM- Newton and Suzaku observations 
have also been used to study the Si/Fe ratio 



(e.g. Matsushita et al.l 120031 iTamura et al.l 


2004 


de Plaa et al.l 120071: |Rasmussen & Ponman 


2007 


Matsushita et al.ll2007a b; ISato et al. 12007, 


2008 


Komiyama et al. 2009; Simionescu ct al. 


2009 


de Grandi fe Molendil 120091: ISato et al.l l2009allb 


2010: Sakuma et al.1 l201lh. With Suzaku obser- 



vations of clusters and groups with ICM temper- 
atures lower than ~ 4 keV, the derived values 
of Si abundance are close to those of Fe out to 
0.2ri8o-0.3ri8o, with a small scatter when using 
the solar abundance table by lLodders |2003) . On 
the basis of the abundance ratios of Si and Fe, the 
contributions from SN la and SNcc have been esti- 
mated. However, excluding cool-core regions, the 
error bars in the Si/Fe ratio of hotter clusters are 
very l arge (|Tamura et al ] l2004h . Matsushita ct al 



(|2012f ) derived the Si/Fe ratio in the ICM of the 
Coma cluster from the flux ratios of the Lya line 
of H-like Si and the Ka line of He-like Fe. The 
temperature dependence of the line ratio above 
several keV is relatively small and limits the sys- 
tematic uncertainty in the derived abundance ra- 



tio. The derived Si/Fe ratio of the Coma cluster is 
approximately 1 solar according to the same solar 
abundance table, with no radial gradient out to 
0.5ri 80 . 

The Perseus cluster is the brightest cool-core 
cluster with a redshift of 0.018. With Chandra 
observations, the cool core of the Perseus clus- 
ter shows complicated features su ch as ripples 



and shocks around radio bubbles (Fabian et al 



20061 ). With XMM, ITamura et ail ((2004J) derived 



the Si/Fe ratio outside the cool core of the Perseus 
cluster, (50-200 h^ kpc), to be 0.77 ± 0.25 in 
solar un i ts usi ng the solar abundance table by 
Lodderd ( 2003 ). Here, the Hubble constant is 
km s _1 Mpc" 1 . With Suzaku data, 



WOh 



100 



Tamura et al. (l2009l) derived the abundance dis- 



tribution of Mg, Si, S, Ar, Ca, Cr, Mn, Fe, and 
Ni of the centra l regio n of the Perseus cluster. 



Simionescu et al.l (|2011l ) derived the electron den- 
sity and Fe abundance profiles out to the virial ra- 
dius toward the northwest (NW ) and east (E) di- 
rections observed with Suzaku. ISimionescu et aT 



(|2012l ) studied the large-scale motions of the ICM 
of the Perseus cluster using ROSAT, XMM, and 
Suzaku data. 

In this paper, we study the Si/Fe and S/Fe ra- 
tios in the ICM of the Perseus cluster observed 
with XMM out to 1.1 Mpc. We collect the K- 
band luminosity of galaxies and calculated the 
IMLR and Si-mass-to-light ratios (SMLR) out to 
this radius. We also calculated the IMLR and con- 
strained SMLR out to 1.8 Mpc, or about the virial 
radius, using the Suza ku measurements in th e 
NW and E directions by ISimionescu et ah ( 2011 ) . 
Using a Hubble constant of 70 km s _1 Mpc -1 , 
1' corresponds to 22 kpc. The value of ri 80 
expected from numerical sim ulations is nan = 
1.95 Knax/k(T)l\Q keV Mpc (jEvrard et alJll99fit 



Markevitch et al.lll998l ). For the Perseus cluster, 
this value corresponds to 2.2 Mpc using the aver- 
age temperature fc(T)=6.1 keV. Hereafter, we de- 
note this radius as ri8o<kT>- An fit with the uni- 
versal m ass model of Navarro , Fren k and White 



(NFW) ({Navarro et al.l Il996l 119971 ) to the hy- 



drostatic mass derived from the Suzaku obser- 
vations in the NW dir ection gave r2nn (herea fter 
^200<he>) as 1.8 Mpc ( Simionescu et al.ll2011 ). 

We used the solar abundance table by lLodders! 
(|2003l ). in which the solar Si, S and Fe abun- 
dances relative to H are 3.47xl0~ 5 , 1.55 xlO -5 , 
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Right ascension 



0.02 0.04 0.08 

Fig. 1. — Exposure-corrected combined MOS im- 
age of the Perseus cluster (0.5-4.0 keV). The cir- 
cles have radii of 0.03, 0.06, 0.1, 0.2, 0.3, and 
0.5ri8o<kT>- The numbers below the color bar 
have units of counts s pixel . 

and 2.95xl0~ 5 , respectively, by number. Consid- 
ering a difference in solar He abund ance, the Fe 
abundance yielded bv lLoddersI (l2003h is 1.5 times 
higher than that us i ng the photospheric value by 
Anders fe Grevessel (|l989t ). Using the table by 
Loddersl ( 2003 ). the Si/Fe and S/Fe ratios are fac- 
tors of 1.5 5 and 1.52, respect i vely, smaller than 
those from Anders fe Grevessel (|l989h . 

In this paper, errors were quoted at a 68% con- 
fidence level for a single parameter. 

2. Observations 

We analyzed the archival data of the XMM- 
Newton observations of the Perseus cluster using 
PN, MOS1, and MOS2 detectors. In this study, 
we used SASvl2.0, but the details of observations, 
event selection, and background subtraction are 
the same as those in Matsushita (l201lh . The ob- 
servation log is summarized in Table [TJ We in- 
cluded three observations in the archive in addi- 



tion to those in lMatsushital ([201 if ). The exposure- 
corrected combined MOS image of the Perseus 
cluster within an energy range of 0.5-4.0 keV is 
shown in Figure [TJ Spectra were accumulated in 



concentric annular regions centered on the X-ray 
peak of the Perseus cluster. The spectra from 
MOS1 and MOS2 were added. 

3. Data Analysis and Results 

3.1. Abundance ratios of Si/Fe and S/Fe 
out to 0.5ri 80 <kT> 

Figure [2] shows representative MOS and PN 
spectra around the Lya lines of H-like Si and S 
and Ka line of He-like Fe. Although the Si and S 
lines are clearly seen in the spectra, the peak level 
of the lines in the spectra is only several percentage 
points above that of the continuum. For the spec- 
tra beyond several hundred kpc, around the Fe 
line, most of the continuum comes from non- X-ray 
background (NXB). As a result, a small system- 
atic uncertainty in the response matrix and back- 
ground can cause a large systematic uncertainty 
in the abundance of these elements. For exam- 
ple, when we fitted MOS and PN spectra at 0.1- 
0.2riso<kT> using an en ergy range of 0.5-1 0.0 keV 



with a vAPEC model ( Smith et al.l 120011 ). there 



are discrepancies of a few percent around the Si 
and S lines between the data and the model and 
the derived abundances from the PN and MOS 
detectors differ by a factor of two. When we re- 
stricted energy ranges around H-like lines of Si and 
S and refitted the spectra, the derived abundances 
strongly depend on the adopted energy range. 

To derive the strength of the Si line accu- 
rately and to obtain suitable statistics, we fitted 
the raw (without background subtraction) MOS 
and PN spectra of each annular region outside 60 
kpc (0.03ri 80 <kT>) out to 1090 kpc (0.5ri 80 <kT>) 
simultaneously with a sum of the vAPEC code 
v2.0.1, two Gaussians for the Ly a line of H-like Si 
and the Ka line of He-like Si, a power-law model 
with r = 1.4 for the cosmic X-ray background 
(CXB), and a power- law model that is not folded 
through the auxiliary response file (ARF) for the 
NXB. Here, we used an energy range of 1.8-2.1 
keV, the Si abundance of the vAPEC model was 
fixed at 0, and the abundance of the other met- 
als was assumed to have a same value which was 
left free. Above 1.8 keV, a systematic uncertainty 
caused by a strong instrumental fluorescence line 
at ~1.7 keV of the MOS detector does not affect 
the derived Si abundance. To derive S and Fe 
line strengths, we fitted the spectra in the same 
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Table 1 

Observation log of the Perseus cluster 



obsid a (RA, Dec) b Date c Exposures (ks) d 



0085110101 


03 h 


19 11 


1 48.36 s 


+41° 


30' 


40.6" 


2001-01-30 


49, 


51, 


48 


0085590201 


03 h 


19 11 


1 45.12 s 


+41° 


05' 


02.5" 


2001-02-10 


42, 


40, 


38 


0151560101 


03 h 


16 11 


1 39.07 s 


+41° 


18' 


38.5" 


2003-02-26 


25, 


25, 


18 


0204720101 


03 h 


21 11 


1 33.80 s 


+41° 


31' 


03.1" 


2004-02-04 


15, 


15, 


12 


0204720201 


03 h 


23 11 


1 18.92 s 


+41° 


30' 


57.2" 


2004-02-04 


25, 


25, 


21 


0305720301 


03 h 


22 11 


1 16.00 s 


+41° 


11' 


28.0" 


2005-08-03 


20, 


21, 


16 


0305720101 


03 h 


18 ff 


1 01.94 s 


+41° 


46' 


46.2" 


2005-09-01 


13, 


13, 


10 


0305690101 


03 h 


1? ir 


' 58.31 s 


+41° 


16' 


15.7" 


2006-02-10 


26, 


27, 


26 


0305690301 


03 h 


19 11 


1 45.91 s 


+41° 


52' 


51.0" 


2006-02-11 


20, 


20, 


16 


0305690401 


03 h 


21 11 


1 49.07 s 


+41° 


48' 


47.2" 


2006-02-11 


27, 


27, 


25 


0405410101 


03 h 


21 11 


1 04.82 s 


+41° 


56' 


46.7" 


2006-08-03 


18, 


16, 


12 


0405410201 


03 h 


18 11 


1 44.75 s 


+41° 


07' 


11.2" 


2006-08-03 


17. 


, 26. 


, 9 



a XMM-Newton observational identifier b in J2000.0 c Date of start of the observation d Exposure times of 

MOS1, MOS2, and PN, respectively, after screenings 



Table 2 

Results of the spectral fitting around Ka line of Fe 



radius 




kT b 


X 2 /d.o..r 


Fe d 


(kpc/ri 80 <A;T>) 




(keV) 




(solar) 


60-130/0.03-0.06 


0.23 ± 0.01 


5.9 ± 0.1 


76/93 


0.63 ± 0.01 


130-220/0.06-0.1 


0.30 ± 0.02 


6.5 ± 0.2 


77/93 


0.54 ± 0.03 


220-440/0.1-0.2 


0.27 ± 0.03 


6.2 ± 0.2 


128/124 


0.47 ± 0.03 


440-650/0.2-0.3 


0.23 ± 0.05 


5.9 ± 0.5 


59/52 


0.38 ± 0.04 


650-1090/0.3-0.5 




5.6 ± 0.6 C 


45/52 


0.48 ± 0.07 



a Ratio of flux in units of photons cm _2 s _1 of Lya line of H-likc Fc and Ka line of He-like Fe. b The ICM 
temperature derived from the Fe line ratios using the theoretical expectation by APEC v2.0.1. c x 2 and 
degrees of freedom of the spectral fitting around the Ka line of He-like Fe, or 6.0-7.2 keV. d Fe abundance 
derived from the flux ratio of Ka line of He-like Fe and the continuum at 3.5-6.0 keV. ° ICM temperature 
derived from the spectral fitting within energy range of 0.8-7.3 keV, considering 10% systematic error. 
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3x10' 




2.3 2.4 2.5 2.6 2.7 2.3 2.4 2.5 2.6 2.7 

Energy (keV) Energy (keV) 




Energy (keV) Energy (keV) 

Fig. 2.— Spectra at 440-650 kpc or 0.2-0.3ri 80 < k T> (left panels) and 650-1090 kpc or 0.3-0.5ri 80<k T> 
(right panels) around Lya lines of H-like Si (top panels) and S (middle panels) and Ka line of He-like Fe 
(bottom panels). Open (black) and closed (red) circles correspond to the MOS and PN spectra, respectively. 
Here, the background was not subtracted. The solid lines correspond to the best-fit model, which is a sum 
of a vAPEC v2.0.1 model, Gaussians, and background components. The dashed lines show the contribution 
of the continuum of the ICM component and the background. The dot-dashed lines in the bottom panels 
show the contribution of the Fe lines from the ICM. The bump at 6.4 keV in the background component 
corresponds to an instrumental line. 



Table 3 

Line ratios of the Ly« of H-like Si and S to Ka to He-like Fe and Si/Fe and S/Fe ratios. 



radius 
(kpc/ri 80 <feT>) 


F Si /F Fc a 


Si/Fe b 
(solar ratio) 


X 2 /d.o.r 


Fs/Fpe a 


S/Fe b 
(solar ratio) 


X 2 /d.o.f d 


60-130/0.03-0.06 

130-220/0.06-0.1 

220-440/0.1-0.2 

440-650/0.2-0.3 

650-1090/0.3-0.5 


0.33 ± 0.03 
0.27 ± 0.04 
0.28 ± 0.03 
0.37 ± 0.05 
0.26 ± 0.08 


0.99 +g£ 

o.87 ±S:ii 

0- 88 ±°fo 

1 - i0 -0.19 
n 77 +0.23 
u - ' ' -0.25 


22/22 
22/22 
16/22 
11/22 
19/22 


0.23 ± 0.02 
0.23 ± 0.03 
0.17 ± 0.03 
0.16 ± 0.05 
0.22 ± 0.06 


1 02 + 009 

1 - uz -0.09 

i.n ±8:iS 
0.76 t° .lt 

0-69 ±gS 
0.97 t° it 


80/71 
49/71 
38/46 
54/46 
45/46 



a Ratio of flux in units of photons cnr 2 s _1 of Lya line of H-like Si or S to Ka line of He-like Fe. b The 

Si/Fe or S/Fe abundance ratios in units of solar ratio derived from the line ratio using the theoretical 
expectation by APEC v2.0.1. c x 2 anCl degrees of freedom of the spectral fitting around the Lya line of 
H-like Si, or 1.8-2.1 keV. d x 2 an d degrees of freedom of the spectral fitting around the Lya line of H-like 



S, or 2.3-2.8 keV. 



manner but used an energy range of 2.3-2.8 keV 
and 6.0-7.2 keV, respectively. We added a 6.4 keV 
Gaussian for an instrumental background line. We 
did not used the region within 60 kpc where the 
temperature structure is rather complic ated be- 
cause of the cool core ( Fabian et al. 20061) . Tables 
[5] and |3] shows the results of the spectral fits. We 
obtained reasonable x 2 values and, as shown in 
Figure [21 the spectra around Si, S and Fe lines 
of the MOS and PN were well reproduced by the 
model. 

Within 0.3riso<kT>, we derived the ICM tem- 
perature from the line ratio of Lya of H-like Fe to 
Ka of He-like Fe using the APEC v2.0.1 plasma 
code to reduce systematic uncertainties caused 
by those in the response matrix and background. 
At 0.3-0.5rigo<kT>, we used deep sky observa- 
tions as a background, and fitted the background- 
subtracted spectra within an energy range of 0.8— 
7.3 keV with an APEC model and a power-law 
without ARF for a possible remaining NXB com- 
ponent. Furthermore, we added an 10% sys- 
tematic error in the derived ICM temperature at 
0.3-0.5ri8o<kT> considering the systematic uncer- 
tainties in the res ponse matrix and background 
(|Matsushital l201ll ). The derived ICM tempera- 
tures are shown in Table [2] 

Figure [3] shows the ratio of Lya of H-like Si 
or S to the Ka line of He-like Fe plotted against 
the ICM temperature. The derived line ratios of 



Si and S to Fe are close to theoretical expecta- 
tions at solar Si/Fe and S/Fe ratios, respectively, 
by APEC v2.0.1 (Figure [3J). These two line ratios 
show a similar temperature dependence: above 5 
keV, the line ratios are relatively flat at fixed Si/Fe 
or S/Fe ratios. The weak temperature dependence 
can minimize the effect of uncertainties in the tem- 
perature structure of the ICM. At a given ICM 
temperature and abundance ratio, the theoretical 
expectations obtained by the MEKAL and APEC 
codes differ by approximately 10 %. The APEC 
vl.3 code gave almost the same line ratios with 
APEC v2.0.1. Therefore, any systematic effect 
due to the plasma codes is expected to be insignif- 
icant. 

Using the APEC v2.0.1 code, we converted the 
derived line ratios to the abundance ratio of Si/Fe 
and S/Fe assuming the single temperature struc- 
ture. Table [3] and Figure 2] summarize the de- 
rived Si/Fe and S/Fe ratios. The derived Si/Fe 
and S/Fe ratios show no radial gradient out to 
~ 1100 kpc, or 0.5riso<kT>. Since the temper- 
ature dependence of the flux ratio of these lines 
is relatively flat above several keV, the derived 
abundance ratios should not change as a result 
of underestimating the ICM temperature. In con- 
trast, if the ICM temperature is overestimated or 
if there is an emission component with a lower 
temperature, the same flux ratio yields a lower 
ICM abundance. As a result, the Si/Fe and S/Fe 
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Fig. 3. — Flux ratios of Ly a line of H-like Si 
(left panel) and S (right panel) to Ka line of He- 
like Fe plotted against plasma temperature. The 
temperatures and line ratios were derived within 
the radial ranges listed in Tables [2] and El Dashed 
lines indicate theoretical ratios with a plasma of 
constant Si/Fe and S/Fe ratios according to APEC 
v2.0.1 model. The numbers on the plots show the 
Si/Fe or S/Fe ratios in solar units. 



ratios should not be much greater than unity in so- 
lar units. The emission- weighted averages of the 
Si/Fe and S/Fe ratios outside the cool core, 130- 
1090 kpc, or 0.06-0.5ri 80<k T> are 0.91 ± 0.08 and 
0.93 ± 0.10, respectively, in units of the solar ratio. 

To derive the Fe abundance, we used the flux ra- 
tio of the Ka line of He-like Fe and the continuum 
at 3.5-6.0 keV (see Matsushita 2011 for details). 
Here, we used the APEC v2.0.1 code and the tem- 
peratures shown in Table [2j The results are shown 
in Table [2] Th e derived Fe abund ances are almost 
the same as in iMatsushital ( 2011 ). Figure [5] shows 
the radial profiles of Si, S, and Fe abundances. Be- 
yond 220 kpc (0.1ri 80 <kT>) from the center, the 
Si and S abundances have flat radial profiles at 
~0.4 solar. 



3.2. if-band luminosity of galaxies 

Because the if-band luminosity of a galaxy cor- 
relates well with the stellar mass, we calculated 
the luminosity profile of the if-band. We collected 
if -band magnitudes of galaxies in a 10 x 10 deg 2 
box centered on the center of the Perseus cluster 
from the Two Micron All Sky Survey (2MASS). 
Figure [6] shows the galaxies detected by 2MASS. 




Fig. 4. — Radial profiles of abundance ratios of 
Si/Fe (closed circles with error bars) and S/Fe (di- 
amonds) in units of solar ratio of the Perseus clus- 
ter, derived from flux ratio of Lya lines of H-like 
Si or S to Ka line of He-like Fe. 



In the Perseus cluster, the distribution of galax- 
ies is elongated in an east-west direction. We 
selected galaxies above the completeness limit of 
2MASS, K s = 13.5 in apparent magnitudcQ-We 
corrected t he foreground Galactic extinction of 
A K = 0.06 (|Schlegel. Finkbeiner. DaTislll998l ) by 
using the NASA/IPAC Extragalactic Database. 
The if-band surface-brightness profile of the se- 
lected galaxies centered on the cD galaxy is shown 
in Figure [7] The average surface brightness in 
the region at 1.0r 180 < kT > < r < 2.0r 180<kT> 
(100' < r < 200') was subtracted as the back- 
ground. Next, we deprojected the brightness pro- 
file assuming a spherical symmetry and derived 
the three-dimensional profile of the if-band lumi- 
nosity. Figure |S] and Table 2] show the integrated 
profile of the if -band luminosity. The total K- 
band luminosity within ri8o<kT> is relatively in- 
sensitive to the background: when we used the re- 
gion 1.0-1.5riso<kT> as the background, the total 
if-band luminosity changed by only 10%. 

The completeness limit of galaxies of 2MASS 
corresponds to the absolute if-band magnitude of 
-21.0. This value is much fainter than the charac- 



1 See http://www.ipac.caltech.edu/2mass/releases/sccond/doc/cxplsup.html 



Table 4 

Cumulative values of ^-band luminosity (L k ), gas mass (M gas ), IMLR, and SMLR derived 

with XMM. 



radius 


L K 




IMLR 


SMLR 


(kpc/ri 80<k T>) 


(1O 12 L ) 


(1O 12 M ) 


(1O- 4 M /L ) 


(10- 4 M Q /i Q ) 


130/0.06 


2.0 


3.2 


12.5 ± 0.2 


7.0 ± 0.1 


220/0.1 


2.9 


6.6 


16.0 ± 0.4 


8.5 ± 0.2 


440/0.2 


4.3 


21 


29 ± 1 


15 ± 1 


650/0.3 


6.0 


39 


34 ± 2 


20 ± 1 


1090/0.5 


6.9 
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Fig. 5. — Radial profiles of abundances of Si 
(closed circles with error bars), S (solid diamonds), 
and Fe (dashed diamonds). 



teristic magnitude, of galaxies in the P erseus 
cluster of -25.09 ± 0.29 (|Lin fe Mohrll2004h . Inte- 
grating the luminosity function with the assump- 
tion of the Shechter function, the contribution of 
fainter galaxies below the 2MASS limit is a few %. 



3.3. Gas-mass, Si-mass, and Fe-mass pro- 
files 

To estimate the integrated gas- and Fe-mass 
profiles, we accumulated MOS spectra in concen- 
tric annular regions with width of l'-2' out to 50', 
or 0.5ri 80 <fcT>- We subtracted those of deep sky 
observations accumulated over the same detector 
regions as background. We fitted these spectra 
within an energy range of 1.6-5.0 keV with an 




Fig. 6. — Distribution of galaxies detected by 
2MASS in the -ftf-band. The circles correspond to 
1.0 and 2.0ri 80 <kT>. 



APEC model to avoid uncertainties in the back- 
ground. Here, the temperature and Fe abundance 
were restricted within the error bars shown in 
Table [5] We fitted the radial profile of the de- 
rived emission measures from the spectral fits per 
area with a sum of two /3-models. The profile 
is well represented with the model and we calcu- 
lated the electron density profile. We integrated 
the electron density profiles observed with XMM 
in this way out to 0.5nso^kT> and Suzaku to- 
ward E and NW directions by ISimionescu et al 
(J201l|) out to ~ T2oo<he> and derived the in- 
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Fig. 7. — Surface brightness profile of galaxies 
around the Perseus cluster observed with 2MASS 
photometric data in the -ftT-band. Dashed lines 
correspond to 1.0riso<kT> and 2.0ri8o<kT>- The 
dotted line shows the adopted background level. 

tegrated gas-mass profiles (Figure [5] and Table 
IU). Because of gas sloshing and cold front at 700 
kpc in the E direction, the integrated gas mass 
toward this direction is higher than that of the 
NW direction. The gas-mass profile with XMM 
falls between the profiles of the two directions 
with Suzaku out to 0.5ri8o<kT>- With the inte- 
grated if-band luminosity and gas mass profiles 
with XMM and weighted average for the two di- 
rections with Suzaku, we calculated the radial pro- 
files of cumulative gas-mass-to-light ratio. In spite 
of differences in the observed azimuthal directions 
and in the analysis methods, as shown in Figure[9j 
the Suzaku and XMM gave similar radial profiles 
for the gas-mass-to-light ratio. The ratio increases 
with radius out to the virial radius. 

We used the Si and Fe abundance profiles from 
XMM-Newton data and derived integrated mass 
profiles of Fe and Si in the ICM. We also inte- 
grated the Fe mass using the weighted average of 
the Fe abundance profiles for t he E and NW di- 



rections observed with Suzaku ( Simionescu et al 



20111 ). The radial profiles of the cumulative IMLR 
and SMLR are shown in Figure [10] and Table |U 
The error bars of the mass-to-light ratio include 
only the errors in abundance. These profiles con- 
tinue to increase with radius out to 1800 kpc, ~ 
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Fig. 8. — Integrated .KT-band luminosity (thin 
solid line) and gas-mass profiles of the Perseus 
cluster derived with XMM-Newton (bold solid 
line), and Suzaku (E:red dotted line, NW: blue 
dashed line), integrated over the electron density 
profiles, assuming spherical symmetry. The two 
vertical lines show riso<kT> and r2oo<HE>- 



^200<HE>- 

4. Discussions 

4.1. The abundance pattern and contribu- 
tion of SN la and SNcc 

The derived Si/Fe and S/Fc ratios in the ICM 
of the Perseus cluster show no radial gradient 
out to 1.1 Mpc, or 0.5riso<kT>- The emission- 
weighted averages of the Si/Fe and S/Fe ratios 
within 0.06-0.5ri 80<k T> are 0.91 ± 0.08 and 0.93 
± 0.10, respectively, in solar units. These values 
are consistent with the emission- weighted average 
of the Si/Fe ratio of 0.99 ± 0.13 in solar units 
of the Coma cluster within 0.5ri8o<kT>- They 
are also consistent with the values of ~ 0.8 in so- 
lar units of those in clusters and groups observed 
with Suzaku with ICM temperature le s s than sev- 



eral keV 



e.g. 



Matsushita ct al. 2007a; Sato et al 



20071 12008L l2009allbl . 1 20 id: iKomivama et al.l l2009; 



Sakuma et aTll201ll ). 

Based on the abundance ratios including Si and 
Fe, the contribution from SN la and SNcc were de- 



rived ( e.g.. iFin ogucno v et al 
20071 : iRasmussen fe Ponmanl 



l200d Ide Plaa et al 



20071 ISato et al 
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Fig. 9. — Radial profiles of integrated gas-mass- 
to-light ratio with XMM (solid line) and Suzaku 
(dashed line). Here, we used the weighted average 
of the E and NW results of Suzaku. 



20071: Ide Grandi fe Molendl2009t lSi mioiiescu etal 



20091: iMatsushita et al.ll2012t ). Since Si and Fe 
were synthesized by both SN la and SNcc, the 
contribution of the two types of SN strongly de- 
pends on the adopted nucleosynthesis model of 
SN la. With the yields of S Ncc w ith metallic- 
ity = 0.02 by iNomoto et ail (|2006h . they found 
that by using the clas sical deflagration model, W7 
(jlwamoto et all I1999IK and a delaye d detonation 
(DD) model, WDD3 jlwamoto et aljl999h . for the 
theoretical SN la yield, over a half of the Fe and a 
few tens of percent of Si in the ICM were synthe- 
sized in SN la. By adopting another DD model, 



ag ; 

WDD1 ()Iwamoto et al.lll999h . they found that ap 



proximately half of the Si and most of the Fe come 
from SN la. Suzaku enabled us to measure O and 
Mg abundances in the ICM outside cool cores of 
clusters and groups of galaxies with ICM temper- 



Sato et all (|2007l 



atures small er than several keV 
2009aL l2010h found that the mixture of yields of 
the W7 model and SNcc gave better fits of the 
observed abundance pattern of O, Mg, Si, S, and 
Fe in the ICM observed with Suzaku than for the 
WDD1 model and SNcc. The number ratio of 
SNcc and SN la to synthesize metals in the ICM 
was a lso estimated with Suzaku and XMM data 
(e.g. Ide Plaa et al 



2007; de Grandi fe Molendi 
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Fig. 10.— Radial profiles of IMLR with XMM 
(black closed circles) and Suzaku (green open cir- 
cles) and SMLR with XMM (red crosses). Here, 
we used the weighted average of the E and NW re- 
sults of Suzaku and assumed spherical symmetry. 
The orange star indicates the SMLR assuming the 
Si/Fe ratio beyond 0.5ri 80 <kT> is 0.91 ± 0.08 in 
solar units. 



2009: ISimionescu et al.1 120091: ISato et all 12007 



l2009al |2010| ). Using the Suzaku results includ- 
ing O and Mg, the number ratio of SNcc to SN la 
was estimated to be ~ 3.5 and ~ 2, using the W7 
and WDD1 models, respectively. 

We also calculated yields mixtures of nucleosyn- 
thesis models. By adopting a Si/Fe ratio of 0.91 ± 
0.08 in solar units in the ICM of the Perseus clus- 
ter, and using the W7 model, we find that 65-74 % 
and 23-30 % of Fe and Si, respectively, were syn- 
thesized by SN la. Using the WDD1 model, 80-91 
% and 56-75 % of Fe and Si, respectively, origi- 
nated from SN la. Here, we also used the yields 
of SN cc with metallicity = 0.02 by INomoto et al. 
(|2006l K The difference in the Si/Fe ratios in the 



yields of SNcc assuming a Salpeter initial mass 
function (IMF) and a top-hea vy IMF and differ 



ent nu cleos ynthesis models by I Chieffi fe Limongi 
(|2004l ) and IWooslev fe Weaver! (|l995h are rekv 
tively small, within ~10%. 

From the Si/Fe ratio of the Perseus cluster, we 
also derived the number ratio of the SNcc to SN 
la to contribute metals in the ICM. Using the W7 
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and WDD1 models, the derived number ratios are 
3.2-4.8 and 0.8-2.0, respectively. The number ra- 
tio from the Si/Fe ratio strongly depends on the 
adopted nucleosynthesis model, as found by pre- 
vious studies. 

The average S/Si ratio of the ICM in the 
Perseus cluster is 1.02 ± 0.14 in solar units. 
Th e S/Fe ratio o f the theoretical yield of SNcc 
by iNomoto et ail (j2006l) with metallicity = 0.02 
and the Salpeter IMF is 0.83 in solar units. The 
yields of the W7 model and the WDD models are 
approximately 1.1-1.2 in solar units. Therefore, 
the observed S/Si ratio of the Perseus cluster is 
consistent with the yields of both SN la and SNcc. 

4.2. Si-Mass-to-Light Ratios and initial 
mass function of stars 

The integrated SMLR using the i^-band at 
0.5riso<kT> of the Perseus cluster is ~0.003 
Mq/Lk,q- When we assume that the Si/Fe ratio 
beyond 0.5riso<kT> is the same as that within 
the radius, the SMLR out to 0.86ri 8 o<kT> be- 
comes 0.004-0.005 M Q /L KiQ (Figure Q33]) . The Si 
abundance is not expected to increase with radius, 
because gas is more extended than stars. Then, 
assuming the Si abundance beyond 0.5rigo is 0.4 
solar, which is the Si abundance at 0.1-0.5riso, the 
SMLR becomes 0.005 M ^JL K ^. If the clumping 



of the gas is significant (jSimionescu et al.l 120111 ) 
and the gas fraction at the virial radius is close to 
the value of the Wi lkinson Microwave An isotropy 
Probe (WMAP) 7 (jKomatsu et al.l 120 111 ), the cu- 
mulative SMLR out to 0.86ri8o<kT> is overesti- 
mated by 20-30%. To summarize, the expected 
value of the SMLR at 0.86 ri 80<k T>, is 0.003- 
0.005 M Q /L K<& . Adopting the W7 and WDD1 
models for the SN la nucleosynthesis model, the 
SMLR synthesized by SNcc becomes 0.002-0.004 
Mq/L K}Q and 0.001-0.002 M e /L KiQ , respec- 
tively. 

The O, Mg, and Si abundances in the hot 
interstellar medium (ISM) in early-type galax- 
ies reflect stellar metallicity because the ISM is 
thought to come from stellar mass loss. The 
abundances of these elements in the hot ISM in 
bright early-type galaxies observe d with Suzaku 
are approximately 0.5-2 solar ([Konami et al. 
20121 ). Extrapolating the observed gradient of 



the mean stellar metallicity, [Fe/H] of individual 
galaxies ranges from -0.8 to +0.3 and correlates 
with stellar luminosity. They also found that typ- 
ical [Mg/Fe] is approximately +0.2. These re- 
sults indicate that the Si abundance in stars in 
giant cluster galaxies is approximately 0.5-2 so- 
lar. Using i-T-band, the stellar mass-to-light ratio 
in bri ght early-type galaxies are approximately 
unity (jNagino fe Matsushita! 120091 ) . As a result, 
the SMLR trapped in stars should be 0.0005- 
0.002 Mq/Lk,q- The estimated value of the total 
SMLR in the Perseus cluster (i.e., the sum of the 
SMLR in the ICM and in stars) is 0.004-0.008 
M q /L k ,q- 

Theoretical models predict that the oxygen- 
mass-to-light ratio (OMLR) and SMLR of a clus- 
ter ar e very sensitiv e functions of the slope of the 
IMF ( Renzini 20051 ). Here, the oxygen and silicon 
mass are a sum of that trapped in stars and that 
in the ICM. Adopting a Salpeter IMF with a slope 
of 2.35, and difference in the stellar mass-to-light 
ratio between the B-band and ii"-band equal 5, we 
find that the expected value of the SMLR from 
SNcc is ~ 0.002 Mq/L k ,q- This value is close 
to the sum of the SMLR in stars and ICM from 
SNcc, adopting the nucleosynthesis yield of WDD1 
model. Using the W7 yields, the SMLR from SNcc 
correspo nds to a slope o f ~ 2 based on the calcu- 
lation bv lRenzinil (|2005l ). A top-heavy IMF with a 
slope of 1.35 overproduces metals more than that 
with a factor of 20. Therefore, the expected value 
of the SMLR out to the virial radius does not need 
the top-heavy IMF. 

4.3. Fe mass and past SN la rate 

The solar Si/Fe abundance ratio in the ICM 
of the Perseus cluster indicates that most of the 
Fe was synthesized by SN la. The estimated rate 
of the current SN la rate in present early-type 
alaxies are 0.1-0.3 SN Ia/(100 yT)/(10 10 L Rr7l ) 
Cappellaro et al.lll997lll999HTuratto et al.l ll999; 



Sharon et al.l l2007t iMannucci et al.l 120081 ). We 

adopted the Fe mas s per SN la rate by W7 



the Mg index in optical spectra o f elliptical galax- 
ies, Kobayashi fe Arimotol ( 19991 ) calculated that 



model of 0.75 Mq (jlwamoto et al.1 1999?) and 
Lk/Lb~ ^>Lkc^I Lbi^ for early- type galaxies 
([Nagino fe Matsushita! l2009h . Then, accumu- 
lating the present SN la rate over the Hubble 
time, 13.7 Gyr, the expected IMLR from the SN 
la becomes (2-6)x1O^ 4 M /_Lk, q . This value 
is over an order of magnitude smaller than the 
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observed IMLR (6-7) xlO- 3 M e /L K<cs and (7 
9)x10~ 3 Mq/L k ,q within 0.5 ri 80 <kT> and ~0.86 
r-i80<kT> , respectively. 

The increase of the radial profile of the IMLR 
with radius of clusters indicates that Fe in the 
ICM extends farther than stars, at least out to 



0.5 



Leccardi fe Molendil (J2008) and 



4>r- 

Matsushital (|2011l ) discovered that the Fe abun- 
dance profiles of the ICM within rsoo are flat- 
ter thaii_exp_ected from the numerical simulations 
bv iFabian et all (|2008l ). without AGN feedback. 
These results indicate that a significant fraction of 
Fe is synthesized in an early phase of cluster evo- 
lution, because if metal enrichment occurs after 
the formation of clusters, the metal distribution is 
expected to follow the stellar distribution. Consid- 
ering that most Fe is synthesized by SN la, the life- 
times of most of SN la are much shorter than the 
Hubble time, and the SN la rate in cluster galaxies 
was much higher in the past. If the IMF is close 
to that in our Galaxy and if most of stars in our 
Galaxy and clusters were already formed before a 
few Gyrs ago, the abundance pattern (including 
the ICM and stars in cluster galaxies) should nat- 
urally be similar to the solar abundance pattern. 

4.4. Comparison of radial profiles of IMLR 
with other systems 

Figure [TT] compares the cumulative IMLR pro- 
file of the Perseus cluster with that of clusters 
of galaxies observed with Suzaku or XMM out 
to 0.5-1.0ri 80 <kT> including Coma (8 keV; Mat- 
sushita et al. 2012), Abell 262 (2 keV; Sato et 
al. 2009b), AWM 7 (3.6 keV; Sato et al. 2008), 
and the Hydra A (3.0 keV; Sato et al. 2012) clus- 
ters. The IMLR profiles of three cool-core clusters 
at a similar redshift, Abell 262, AWM 7, and the 
Perseus cluster agree very well with each other at 
0.2-0.5 riso<kT>, whereas the Coma cluster has 
a lower IMLR at 0.5 ri 80 <kT>- There is no sig- 
nificant temperature dependence on the derived 
IMLR profiles. Since the IMLR profiles increase 
with radius, comparison of the total IMLR values 
requires observations of IMLR profiles of clusters 
out to the virial radii. 

At a given radius in units of ri8o<kT> , the 
Hydra A cluster has a factor-of-two-higher IMLR 
than the Perseus cluster. Because of the higher 
redshift of the Hydra A cluster and lower ICM 
temperature, the number of galaxies detected in 
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Fig. 11. — Profiles of IMLR plotted against ra- 
dius scaled with ri 8 o<kT> of the Perseus clus- 
ter (closed squares with XMM and op en squares 
with S uzaku), AWM 7 (open triangles. ISato et aL 
(|2008l ), Sato et al. in preparation) , Abell 262 



(open diamonds, Sat o et al 



(2009t])j, Hydra A 



duster (crosses, E^IIIiS&th; Coma 



cluster (closed circles, Matsushita et al.l ( 20121) ). 



Hydra A might not be sufficient to deprojected 
the luminosity profiles in K-hanA. However, at 
~ r i80<kT> ! where systematic uncertainties due 
to the limited number of galaxies are relatively 
small, the cumulative IMLR of the Hydra A and 
the Perseus clusters agree well. No systematic de- 
pendence on the ICM temperature is evident for 
these two clusters with temperatures of 3 keV and 
6 keV. 

At the center, within 0.1 ri8o<kT>, the Perseus 
cluster and the Hydra A have much higher IMLR 
than the other systems, which is due to the very 
bright cool core of these two clusters. Within this 
radius, the cD galaxies dominate the fC-band lu- 
minosity and the systematic uncertainties in the 
K-band luminosity should be small. The IMLR 



12 



within 0.1 rigo<kT> of the Perseus and Hydra A 
clusters are about 1.5-2 xlO -3 Mq/Lk,q, which 
is a factor 3-10 higher than the expected IMLR 
from SN la assuming the present SN la rate over 
the Hubble time. The higher IMLR in these 
two clusters reflect higher gas-mass-to-light ra- 
tio in 0.1 riso<kT>, since the Fe abundance of 
these two clusters are not hig her than those of 



the other clusters in th is region ( Simion escu et al 
20091: lMatsushitall201lh . ISimionescu et al.l (|2009h 



measured the distribution of O, Si, S, Ar, Ca, Fe 
and Ni in the ICM of the Hydra A cluster within 
~0.1 ri8o<kT>- They found O abundance de- 
creases with radius and amount of metals are much 
higher than the expected values by stellar winds. 
They discussed the initial enrichment by SNcc in 
the early phase in cluster evolution and mixing of 
these metals. The O/Fe and Mg/Fe ratios within 
the cool core of the Perseus cluster are about 



unity in solar units (jMatsushita fe Tamura 11201 lh 
and agrees with the O/ Fe ratio in the cool cor e 
of the Hydra A cluster ( Simionescu et al. 2009). 
The metals already synthesized in the protoclus- 
ter phase can dominate those in the central regions 
of these two clusters. 

4.5. Gas-mass-to-light ratio 

The integrated gas mass to light ratio of 
the Perseus cluster increases with radius out 
to 1800 kpc, or ~ r 2 oo<HE>- At 1800 kpc, 
M gas (< r)/L K {< r) = 16 M Q /L Q . This value is 
not so different from 21 ± 4 M Q / L Q for the Hydra 
A cluster (< kT >—3 keV) at r 2 oo<HE>- The stel- 
lar and gas mass fract ions within rsnn depend on 
the total system mass dLin. Mohr. Stanford Ife 003; 



Lin fe Mohj |2004 ISun et all 120091 : iGiodini et al . 
20091 : IZhang et alll2010f ). The gas density profiles 



in the central regions of groups and poor clusters 
were observed to be shallower than those obtained 
by the self-similar model, and the relative en- 
tropy level was correspondingly higher than that 
in rich clusters dPonman. Cannon. Navarro" 19991: 



Pon man. Sanderson. Finoguenovl 12003c iSun et aL. 
2009). Then, the difference in the ratio of gas 
mass to stellar mass might reflect differences in 
distributions of gas and stars, which reflects the 
history of energy injection from galaxies into the 
ICM. To study the fractions of stars and gas in 
clusters of galaxies, and their dependence on the 
system mass, measurements beyond rsoo of other 



clusters are required. 

Based on the Suzaku observations of the 
NW and E dire c tions of the Perseus cluster, 
Simionescu et all ([201 ll ) proposed that the gas- 
clumping effect is significant beyond t^qq, and 
that the electron density is overestimated. As a 
result, Mg as within T2oo<he> is overestimated by 
a factor of 1.3-1.5. Adopting this value for the 
gas mass, the gas-mass-to-light ratio at r2oo<HE> 
becomes 11-12 Mq/Lq, which is close to the value 
at 1000 kpc, or 0.6 T2oo<he>- In other words, if 
there is a significant clumping effect, the gas-mass- 
to-light ratio becomes fiat beyond this radius out 

to r 2 00<HE>- 

5. Summary and Conclusion 

We analyzed XMM-Newton data of the Perseus 
cluster out to 0.5rigo<kT> and derived the Si/Fe 
and S/Fe ratios in the ICM from the flux ratios 
of the Ly a lines of H-like Si and S to the Ka 
line of He-like Fe. The small temperature depen- 
dence of the line ratio limits the systematic uncer- 
tainty in the derived abundance ratio. The derived 
Si/Fe and S/Fe ratios in the ICM show no radial 
gradient. The emission- weighted averages of the 
Si/Fe and S/Fe ratios beyond the cool core of the 
Perseus cluster are 0.91 ± 0.08 and 0.93 ± 0.10, 
respectively, in solar units. These abundance ra- 
tios indicate that most of Fe in the ICM within 
0.5riso<kT> was synthesized by SN la. 

We collected if-band luminosities in galaxies 
detected with 2MASS and derived the cumulative 
radial profile of the if -band luminosity of stars 
in the Perseus cluster. We calculated the cumu- 
lative IMLR and SMLR profiles out to ri 8 o<kT> ■ 
Furthermore, by using the electron density profiles 
and Fe abundance profiles in the two directions 
observed with Suzaku, we calculated the cumula- 
tive IMLR profile out to the virial radius. We also 
constrained the SMLR value out to this radius. 
The SMLR of the Perseus cluster out to the virial 
radius is significantly smaller than expected from 
the top-heavy IMF, which has a slope of 1.35, and 
is more consistent of a slope of the Salpeter IMF. 
With the present SN la rate, the IMLR within the 
virial radius is an order-of-magnitude higher than 
the SN la yield accumulated over the Hubble time. 
Since the IMLR increases with radius, the most of 
Fe from SN la were synthesized in an early-phase 
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of cluster formation. 
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